Electron beam formation from spin-orbit interactions in zincblende semiconductor 

quantum wells 
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We find a dramatic enhancement of electron propagation along a narrow range of real-space angles 
from an isotropic source in a two-dimensional quantum well made from a zincblende semiconductor. 
This "electron beam" formation is caused by the interplay between spin-orbit interaction originating 
from a perpendicular electric field to the quantum well and the intrinsic spin-orbit field of the 
zincblende crystal lattice in a quantum well, in situations where the two fields are different in 
strength but of the same order of magnitude. Beam formation is associated with caustics and can 
be described semi-classically using a stationary phase analysis. 



Spin transport in semiconductors [lj-Q can be dramati- 
cally modified by spin-orbit interactions, producing such 
effects as coherent precession without a magnetic fieldp- 
Hoj ]. These effects vanish, however, for transport of unpo- 
larized electrons, as the features generated for an initial 
electron spin polarized up are complementary to those 
generated from an initial electron with spin down. Spin- 
orbit effects have also been found to generate sponta- 
neous spin polarizationfll) and spin currents [l2l - |l6l ] in 
the presence of the flow of unpolarized electrons. An 
open question is whether the spin-averaged transport of 
unpolarized spins can be influenced by the spin-orbit in- 
teraction in a non-trivial fashion. The outward flux 
of electrons injected at a point in a quantum well is 
isotropic (albeit with interesting spin structure), in the 
presence of the spin-orbit field from a perpendicular 
electric field [13, El (Rashba effect [H, HI) or from the 
zincblende lattice of the constituent semiconductors [ilj] 
(Dresselhaus field @, Hi]), or when the two spin-orbit 
fields are of the same strengthfl [Hj. Calculations of 
transport in the mixed Rashba and Dresselhaus system 
based the Boltzmann equation[2^, [24| find isotropic con- 
ductivity for nonmagnetic systems, however these neglect 
interference between states of different momenta. 

Here we describe a dramatic enhancement of electronic 
propagation along a narrow range of real-space angles 
which occurs in the presence of Rashba and Dresselhaus 
spin-orbit fields of specific, different strengths. The an- 
gular width of this "electron beam" depends sensitively 
on the ratio of the strengths of the Rashba and Dres- 
selhaus fields, and the direction of the beam changes 
by 90° when the relative sign of the fields changes. 
This surprising spatial anisotropy, originating from the 
anisotropic dispersion relations of electrons in the two 
fields, is due to general features of the energy con- 
tour surface of the electrons. Furthermore, the elec- 
tron beam formation can be traced, using a station- 
ary phase analysis of the real-space Green's function, 
to coalescing saddle points. Such beams should ap- 
pear in two-contact transconductancefHj as well as other 
transport 26, 27] and scattering [28j] phenomena. A rem- 
nant of this anisotropy (although much weaker) appears 



to cause anisotropy in the dispersion relation of the 
spin[2{| and charge [3fjj| plasmon spectra for a quantum 
well with both Rashba and Dresselhaus fields. 

The conduction band Hamiltonian describing an elec- 
tron confined to a quantum well with [001] growth 
direction [31. 32], to linear order in electron crystal mo- 
mentum, is 

p 2 a 8 
H = 2m + fr^Py ~ a vPx) + -fr(< J xPx - <TyPy)> (!) 

where the first term is the kinetic energy, the second 
the Rashba interaction, and the third is the Dresselhaus 
interaction with strength B. We reparameterize using 



+ /3 2 , /? + 



■exp(zr), (2) 



yielding H = (p 2 /2m) + (Hko/m)U(p), where 

U(p) = sm(T)[a x p y - o-yp x ] + cos(r)[a x p x - a y p y ]. (3) 

U 2 is proportional to the identity matrix: 

U 2 (p) =p 2 (l + S m(2r)sm(2e p ))=p 2 f T (9 p ), (4) 

where p is the angle that the momentum p makes with 
the positive ir-axis. Although U 2 is diagonal in momen- 
tum and independent of spin, it is not isotropic. 

To calculate the electronic propagation we want the 
retarded Green's function in coordinate space, 



G CT>CT ,(r,r') = (r,a| 



1 



E-H 



(5) 



where e — > + . We obtain this from the Fourier transform 
of the momentum space Green's function g CT ,cr'(p, E), 



G a ^{r, r') 
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9a^(p;E). (6) 



Because U 2 is proportional to the 2x2 identity matrix, 
the inverse of E + ie — H can be found directly: 



(2m/h 2 )[(k 2 E -q 2 ) + 2k Q qU(q)} 



2q 2 (k 2 E 



2^o/t) + k E 



(7) 
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FIG. 1: Polar plots of q+(8 q ,T) and q-{6 q ,r), all momenta 
in units of ko- (a) g+, with r = 0.971 and k E = 0.5 (solid), 
1.2 (dashed), 2. (dot-dashed), 4. (dotted), (b) Same as (a) 
but plotting (c) q+, k E = 1.2, r = 0.78 (solid), 0.971 
(dashed), 1.15 (dot-dashed), 1.42 (dotted), (d) Same as (c) 
but plotting q- . Non-convexity for q+ here causes the non- 
convexity in Fig. 4 and produces multiple stationary points. 



where k E = {2mE /h 2 ) 1 ' 2 . 

The zeros of the denominator of the right-hand side of 
Eq.0 q ± = \Q± fc /r /2 |, with Q = (k% + k 2 f T f' 2 , rep- 
resent the intersection of the two energy surfaces E±(q) 
with a constant energy plane, E = (h 2 k 2 E /2m). These 
intersections are shown in Fig. [TJ As the ratio of Rashba 
to Dresselhaus strengths changes (through r) or the elec- 
tron energy changes (through k E ), the constant energy 
surfaces change from convex to non-convex, dramatically 
changing the angular dependence of the electron propa- 
gation. 

Rewriting Eq. ([7]) using partial fractions, 



Gi,i(r) 
2m 



G 2 ,i(r) 



r+T/2 



dd i— [q+hAq+p) + q-h.i(q-p)} (8) 
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i cos(6 q 
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[i+hAq+p) - q-h,i{q-pj\ 



(9) 



where p = rcos(6 q — 6 r ). In these equations the h.j{z) 
are the radial integrals (J pdp) in the Fourier transform. 




2 for various r with positions in units of k 1 
In the upper left panel, r « 7r/4 (correspond- 



FIG. 2: Scaled probability density for unpolarized spins 
rE^|G W (r)," 
and kE ~ 1.2fco. 
ing to a = (3) and the actual variation is less than 1%. In all 
panels the source point is (0,0) (bottom, center). The white 
lines mark the angle where two stationary points coalesce as 
discussed in the text. 



Analytic expressions can be obtained for hj(z) (but are 
not shown as only the asymptotic form J; >m (z) oc exp(zz) 
for large z is relevant here). 

To highlight the angular anisotropy we plot the ra- 
dial distance times the probability density for unpolar- 
ized spins, rY^lj \Gi,j{r)\ 2 , for k E = l-2k Q , in Fig [2 
Only the first two quadrants are shown since the func- 
tion is 7r-periodic. When r = tt/4, corresponding to \a\ = 
|/3 1, there is no anisotropy. As r is increased towards 
Tcritiks), which is defined below (r c rtt (l-2fco) = 1-31169), 
Fig. [2] shows the electron probability pattern narrowing, 
rising and developing considerable structure, finally re- 
solving into a narrow electron beam at t = T cr ;t- When 
t — .871, for example, there are interference patterns 
both in range and angle. The white lines mark a bound- 
ary between the variable, higher probability region and 
a more homogeneous region. The direction of these lines 
is computed below using stationary phase considerations. 
When r w T cr u, the high probability region becomes very 
narrow and intense. Furthermore, along the direction 
9 r = 37r/4, rJ2if \Gij(r)\ 2 increases nearly monotoni- 
cally. Finally when r = —T cr i t , which corresponds to 
changing the sign of the Rashba parameter a while leav- 
ing the Dresselhaus parameter /3 unchanged, the beam is 
reoriented along the direction 8 r — tt/4. The angularly- 
integrated flux for all panels has been confirmed to be 
independent of r. 

Fig. [3] shows the position-dependent polarization, 



P = 



|G M (r)| 2 -|G 2 , 1 (r)| 2 
\G hl (r)\ 2 + \G 2 Ar)\ 2 ' 



(10) 



of an injected spin polarized perpendicular to the quan- 
tum well plane, for r = tt/4, .871, Tcriti and —T cr it when 
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FIG. 3: Polarization, Eq. (|10[) . with positions in units of k 
and k E = l.2k . r = (a) tt/4, (b) 1.15, (c) 1.31 and (d) -1.31. 

ks = 1.2fco- The polarization changes dramatically as 
r changes. The polarization for t = ir/4 suggests the 
spin helix described in Ref. [§]. The shift property 
25- (q + Q) = 2? + (q), used in Ref. [§], yields an ana- 
lytic result for the Green's functions in agreement with 
our results. 

A stationary phase analysis of the angular integrals in 
Eqs. [5] and O provides (1) an explanation for why there is 
a region of high electron probability around the direction 
9 r = 37r/4, (2) the angular width of these regions, and (3) 
the radial dependence of the electron probability along 
the direction 37r/4 for r = T ar n. It also provides a simple 
means of computing the electron probability at large r 
when r and 9 r are not too close to the white lines shown 
in Fig. H 

The key result of stationary phase analysis is that the 
results seen in Fig.[2]result from interfering contributions 
of stationary points, which can coalesce at critical values 
of the spin-orbit field or energy (through r and Ue)- A 
typical integral in the computation is 

= / d0 q %I 1A (q+p) (11) 

where p = rcos(9 q — 9 r ). When r is large, the argument 
of ii.i will vary rapidly with the integration variable, 9 q , 
and dominant contributions to the integral will originate 
at points of stationary phase, where the derivative of the 
argument of I\ A with respect to 9 q vanishes. Apart from 
the endpoints, the argument of I\ A is positive throughout 
the chosen range of integration. Following Ref. [33j, I\ A 
can be replaced by its large-argument asymptotic form 
cx exp(iz). With this replacement, M\ A becomes 

M M ~ f d0 v~zQ~ ex P[*2+ rcos (0g ~ Or)]- 

Figure @] shows polar plots of $ + = <7 + cos(6* — 9 r ). De- 
pending on r and 9 r the phase $ + can exhibit either 
one or three points of stationary phase. The arrows in 
Fig. 3] are drawn along radials which are perpendicular 
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FIG. 4: Polar plots of the phase $+(0 q , r, 6 r ) = 
q+(8 q , t) cos(9q — 9 r ) with 6 r fixed and r fixed. The right hand 
panels show two stationary points coalescing, for r = 0.871 
at r about 35° from 37r/4 and for r = 1.15 at r about 10° 
from 3-7T-/4. The arrows are drawn as normals to the curves 
and correspond the points of stationary phase. Similar plots 
of $_ (not shown) are convex and have only one stationary 
point. 



to the curves, which correspond to points of stationary 
phase. Comparing the polar plots of Fig. [TJ and Fig. 2] it 
is apparent that multiple stationary points arise from the 
non-convexity of the g+ curve. In contrast the g_ curve 
is always convex and therefore the phase $(q~) always 
shows only a single stationary point. The right panels in 
Fig- HI correspond to observation angles (9 r ) such that two 
stationary points are about to merge and disappear leav- 
ing only one stationary point. These observation angles 
correspond to the white lines in Fig. [2] and give a rough 
estimate of the size of the intense areas there. All three 
stationary phase points merge at a value of t — T cr it- 
The lower left panel in Fig. U shows $+ for r near T cr u, 
corresponding to the narrowest beams in Fig. [2] 

At t = T cr it the phase <I> + cx (9 q — 37r/4) 4 . A simple 
stationary phase analysis fails, but the quartic behavior 
implies that M\ t \ cx r^ 1 / 4 , which is consistent with the 
electron probability times r increasing as r 1 / 2 . A radial 
plot (not shown here) of the probability of Fig. [5] along 
37r/4 shows this behavior. Thus the stationary phase 
analysis accounts for the angular widths of the intense 
regions of Fig. [2] they are determined by the transition 
from three points to one point of stationary phase. In 
turn this is where both the first and second derivatives 
of the phase vanish, which is the criterion for the white 
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lines in Fig. [2j We have verified that the angular inte- 
gral of the radial component of the flux times r is in fact 
independent of r. A plot (in EPAPS) of the radial com- 
ponent of the flux times r as a function of 9 r shows that 
the peak flux times r rises proportionally to r 1 / 2 while at 
the same time the peak narrows, as it must for the total 
integral to be fixed. The case of t — it/ 4, corresponding 
to a = f3 is a singular case for which the stationary phase 
arguments here do not apply. 

In the smooth regions of Fig [21 it can be shown that 
there is only a single point of stationary phase associated 
with $ + . This must be combined with the single point 
of stationary phase associated with and the interfer- 
ence of these terms accounts for spin precession and the 
periodicity visible in the polarization plots of Fig. [2] just 
beyond the regions of enhanced probability. 

We conclude with a discussion of length scales for the 
phenomenon of "electron beams" and the implications for 
other measurements in quantum wells. In all our results, 
r, has been shown in units of 1/fco- Experimental values 
of a range from ~ 0.1 — 0.3 eVA in InGaAs quantum 



offers the intriguing possibility of altering and even redi- 
recting the narrow electron beams of Fig. [2] by 90° by 
reversing the sign of the Rashba field (a), by varying an 
electric field perpendicular to the quantum well. 

We acknowledge support from NRI, an ONR MURI 
and helpful discussions with B. Moehlmann. 
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and ^0.5 eVA in In As quantum wells |36|. 
with (3 of the same order. Using the appropriate effec- 
tive masses (~ 0.05m o for InGaAs, ~ 0.025m o for InAs, 
where m is the free mass) this produces a length scale 
of 1/fco ~ 40nm. Thus the scale of the plots, Figs. [2] 
and[3j is about 2/zm. It may be possible to image these 
patterns, as in Ref. [37j . However, the dramatic enhance- 
ments, shown over a longer length scale for clarity, also 
occur on smaller scales ~ 100 nm. For GaAs/AlGaAs a 
is an order of magnitude smaller, meaning the scale of 
Figs. [2] and [3] would be 20/im, however in this system the 
mean free paths frequently exceed 100/xm 38]. 

Fig. [5Js results could be measured directly by 
two-contact transconductance[25j], however the highly 
anisotropic electron propagation dramatically affects 
many other properties. The local density of states near 
impurities has spatial structure described by the square 
of the position-dependent Green's function |28|], and thus 
the wave functions of impurities seen in scanning tun- 
neling microscopy should be highly anisotropic. Scat- 
tering of electrons from impurities will also be highly 
anisotropic, yielding correspondingly anisotropic diffu- 
sive transport [in a real-space Kubo formalism for dif- 
fusive transport [2^, 27 1, the real-space Green's func- 
tions appear in similar combination to Fig. [2] with a 
spatial derivative selecting the conductivity direction 
(a a/ 3 oc J dr[V Q G(r)][V 1 gG(r)])]. These derivatives of 
the Green's function have the same spatial structure as 
shown in Fig. [5J 

We have shown that highly anisotropic electron propa- 
gation (electron beams) occurs in a semiconductor quan- 
tum well at appropriate relative values of the Dresselhaus 
and Rashba spin-orbit fields. In addition to the funda- 
mental consequences for properties of the quantum well, 
driven by the peculiar electronic structure, this system 
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